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Abstract. We assess the cosmological variability of the fine-structure constant a from the analysis of an ensemble of 
FeiU1608, ^2344, ^2374, /12383, A25S1, and /12600 absorption lines at the redshift z = 1.15 toward the QSO HE 0515^414 by 
means of the standard many-multiplet (MM) technique and its revision based on linear regression (RMM). This is the first time 
the MM technique is applied to exceptional high-resolution and high signal-to-noise QSO spectra recorded with the UV- Visual 
Echelle Spectrograph (UVES) at the ESO Very Large Telescope (VLT). Our analysis results in {Aq'/q'>mm = (0.1 ± 1.7) • 10"* 
and (Aa/a)RUM = (-0.4 ± 1.9 ± 2J,y,) ■ 10"*, which are the most stringent bounds hitherto infered from an individual QSO 
absorption system. Our results support the null hypothesis (Aa/a) = at a significance level of 91 percent, whereas the support 
for the result (Aa/a) = -5.7 • 10"* presented in former MM studies is 12 percent. 

Key words, cosmology: observations - quasars; absorption lines - quasars: individual: HE 0515-4414 



1. Introduction 

Modern 10 m class telescopes equipped with instruments 
like the High-Resolution Echelle Spectrograph (HIRES) at 
the Keck Observatory or the UV- Visual Echelle Spectrograph 
(UVES) at the ESO Very Large Telescope (VLT) facilitate the 
accurate observation of QSO absorption (or emission) lines in 
order to study the hypothetical variability of fundamental phys- 
ical constants like the fine-structure constant a = e~ / (AneoHc) 
or the proton-electron mass ratio fi = mp/mg. The interest in 
these studies is motivated by the unification theories incorpo- 
rating vary ing fundamental constants (for a review see, e.g., 
IUzanl2003h . 

From the astronomical point of view the cosmological vari- 
ability of the fine-structure constant is assessed as 



Aa/a - (a; - ao)/ao. 



(1) 



where ao and a- denote the values of the fine-structure con- 
stant in the laboratory and the specific absorption (or emis- 
sion) line system at redshift z, respectively. While observa- 
tional studies based on the fine-structure split ting of intergalac- 
tic aUcali-doublet (AD ) absorption lines (e.g.. lLevshakovll994t 
iMurphv et al.'2001b") and intrinsic QSO emission line doublets 

tBahcall et al 2003.) h ave pro vided robust upper bounds on 

\Aa/a\. iMurphv et a il il2003al hereafter MFW) have recently 
detected a non-zero expectation value (Aa/a) for a sample of 

Send offprint requests to: R. Quast, 
e-mail: rquast@hs . uni-hamburg . de 

* Based on observations made with ESO Telescopes at the La Silla 
or Paranal Observatories under programme ID 066.A-0212. 



143 complex metal absorption systems identified in QSO spec- 
tra recorded with the Keck/HIRES: (Aa/a) = (-5.7+1. 1)-10"^ 
in the redshift range 0.2 < z < 4.2. This remarkable statisti- 
cal evidence for a cosmological variation of the fine-structure 
constant is achieved by means of the many-multiplet (MM) 
technique, which is a generalization of the AD method incor- 
porating the multi-component profile decomposition of many 
transitions from dif ferent multiplets of d ifferent ionic species 
dDzuba et al.Hl999t iMurphv et alJl2003b[ hereafter MWF, and 
references therein). While the MM technique considerably im- 
proves t he formal accuracy , it also shows some immanent defi- 
ciencies. lLevshakovl ( l2003l) illustrates in detail that the prereq- 
uisite assumption that the spatial distribution is the same for all 
ionic species is not valid in typical intergalactic absorption sys- 
tems. Consequently, it is more reliable to apply the MM tech- 
nique to samples of absorption lines arising from only one ionic 
species. Furthermore, in comparison to the AD calculations 
the incorporation of more transitions over a wider wavelength 
range results in a stronger susceptibility to systematic effects. 
Nevertheless, sources of error like wavelength miscalibration, 
spectrograph temperature variations, atmospheric dispersion, 
and isotopic or hyperfine-structure effects do de monstrably not 
explain the detected non-zero expectation value dMurphv et alJ 
1200 lal MWF, MFW). The observational discrepancy grows 
since the radioactive decay rates of certain long-lived nuclei 
deduced from geophysical and meteoritic data provide a strin- 
gent bound, (Aa/a) - (8 + 8) ■ 10"^, back to th e epoch of Solar 
system formation, z < 0.45 dOlive et al.l2003h . 

In the optical spectroscopy, the most striking deficiency 
inherent to all decomposition techniques are non-linear inter- 
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Table 1. Atomic data of the Feii transitions between /1 1608 
and /12600. The laboratory wavelengths A, oscillator strengths 
f, and relativistic c orrection coefficien ts q are excerpted from 
IMWFI Table 2 and lOzuba etaP ilOOt Table 1). The sensitiv- 
ity coefficient Q is defined in Sect. 13.21 Estimiated errors are 
indicated in parentheses 



Tr. 


^(A) 


/ 


q (cm ') 


Q 


1608 


1608.45080(8) 


0.058 


-1300(300) 


-0.021 (5) 


2344 


2344.2130(1) 


0.114 


1210(150) 


0.028 (4) 


2374 


2374.4603 (1) 


0.0313 


1590(150) 


0.038 (4) 


2383 


2382.7642(1) 


0.320 


1460(150) 


0.035 (4) 


2587 


2586.6496(1) 


0.06918 


1490(150) 


0.039 (4) 


2600 


2600.1725(1) 


0.23878 


1330(150) 


0.035 (4) 



parameter correlations preventing the accurate optimization of 
the model parameters and possibly causing ambigous results. 
In fact, the spectral resolution attained in QSO observations at 
the 10 m class telescopes is still not sufficient to resolve the 
metal lines with an expected minimum thermal width of about 
1 kms Clearly, in order to solve the profile decomposition 
problem spectral observations with the highest possible resolu- 
tion and the highest possible signal-to-noise ratio are desirable. 

In this study, we present exceptional high-resolution and 
high signal-to-noise spectra of the notably bright in termediate 
redshift QSO HE 0515-4414 jReimers et alJll998[ z = 1.73, 
B - 15.0) recorded with the VLT/UVES. The spectra reveal 
a multi-component complex of metal absorption lines associ- 
ated with a sub-damped Lyman-cy (sub-DLA) system at the 
redshift z =1-15 (d e la Varga et al. 2000; Ouast et al. 2002^ 
iReimers et al ]l2003l) T We analyze a homogenous subsample of 
Fell .11608, .12344, .12374, .12383, .12587, and .12600 absorp- 
tion lines by means of the standard and a revised MM tech- 
nique in order to assess the cosmological variability of the fine- 
structure constant. 



2. Observations 

HE 0515^414 was observed with UVES during ten nights 
between October 7, 2000 and January 3, 2001. Thirteen ex- 
posures were made in the dichroic mode using standard set- 
tings for the central wavelenghts of 3460/4370 A in the blue, 
and 5800/8600 A in the red. The CCDs were read out in fast 
mode without binning. Individual exposure times were 3600 
and 4500 s, under photometric to clear sky and seeing con- 
ditions ranging from 0.47 to 0.70 arcsec. The slit width was 
0.8 arcsec providing a spectral resolution of about 55 000 in 
the blue and slightly less in the red. The ThAr lamp exposures 
taken immediately after each science exposure provide an ac- 
curate calibration in wavelength. The standard deviation of the 
wavelength versus pixel dispersion solution is about 2.0 mA 
(2.5 mA) in the blue (red), resulting in an absolute accuracy of 
about 0. 15 km s"' in radial velocity space. The raw data frames 
were reduced at the ESO Quality Control Garching using the 
UVES pipeline Data Reduction Software. The calibrate d spec- 
tra we re converted to vacuum wavelengths according to 'Edlen| 
ill 9661) while the barycentric velocity correction was manually 
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Fig. 1. Multi-component Fe ii absorption complex at redshift 
z - 1.15. For convenience, only the transitions .11608 and 
.12383 are shown. Individual components are marked by a ver- 
tical line. The zero point of the radial velocity corresponds 
to the redhsift z = 1.1508. A close-up ranging from -20 to 
100 km s"' is provided in Fig.|2l 

cross-checked using the ESO-Munich Image Data Analysis 
Software and the Image Reduction and Analysis Facility. The 
individual vacuum-barycentric corrected spectra were manu- 
ally cleaned from cosmic rays or pixel defects, rescaled to a 
common median flux level, and resampled to an equidistant 
wavelength grid using natural cubic spline interpolation. The 
combined spectra (Fig.Q show an effective signal-to-noise ra- 
tio per pixel typically better than 100. 

3. Analysis 

3.1. The standard MM technique 

The observed spectral flux F{A) is modelled as the prod- 
uct of the background continuum C(A) and the absorption 
term convoluted with the instrumental profile, i.e., F(A) - 
C(A) J P{^) e"'"(''"f ' d^. While the background continuum is lo- 
cally approximated by an optimized linear combination of 
Legendre polynomials of up to second order, the instrumen- 
tal profile is modelled by a normalized Gaussian given by 
the the spectral resolution of the instrument. Assuming pure 
Doppler broadening, the optical depth t(A) is a superposition 
of Gaussian functions 



NifiAi 



4eomc 



exp 



A-ZjAj 
ZjAibi 



(2) 
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Fig. 2. Close-up of Fig.nshowing the subset of line profiles associated with the central region of the sub-DLA system. The solid 
and dashed curves represent the optimized model and its deconvolution, respectively. The standard deviation of the normalized 
flux is indicated by vertical bars whereas the effective signal-to-noise ratio per pixel is stated explicitly. Data marked by empty 
circles are ignored in the optimization 



where Z; = (1 -h z,)(1 + Uj/c) while Aj, fi, Zi, u,, bi, and A^,- 
denote, respectively, the systemic rest wavelength, the oscilla- 
tor strength, the cosmological redshift, the radial velocity, the 
line broadening velocity, and the column density correspond- 
ing to the line ;. The systemic rest wavenumber w, = is 
parametrized as 

Ui - Wo,/ + qi (Aa/a) (2 + Aa/a), (3) 

where cjoj - l/zlo,, is the wavenumber in the laboratory a nd 
is the relativistic coiTection coefficient jPzuba e t al.'2002'). 

Even though the sub-DLA system exhibits many additional 
metal absorption lines (Al ii, Al iii, Mg i, Mg ii. Si ii, Cr ii, Ni ii, 
and Znii) typically incorporated in the standard MM analysis, 
we decide to apply this technique to the Fe ii transitions only. 
The restriction to one ionic species avoids systematic effects 
if the spatial distribution is not the same for all species, and 
reduces systematic effects arising from isotopic line shifts. In 
addition, the set of Pen transitions between /1 1608 and /12600 
already provides a very sensitive combination for probing the 
variability of the fine-structure constant. 

Each absorption component is modelled by a superposition 
of Doppler profiles with identical radial velocities, widths, and 
column densities. In addition, the value of Aa/a is confined 
to be the same for all components. In order to find the opti- 
mal set of parameter values for both t{A) as well as C{A), the 
weighted sum of residual squares is minimized by means of an 



evolution stra tegy (ES) based on t he concept of covariance ma- 
trix adaption ( iHansen et al.ll2003l) . We ignore the ensemble of 
Fell /12383 and /12600 lines associated with the central region 
of the sub-DLA system since the profiles are saturated and are 
otherwise overemphasized in the optimization. 

3.2. The regression MM (RMM) technique 

The standard MM technique can conveniently be revised to 
avoid the deficiencies pointed out in Sect. \l\ In fact, this 
technique is essentially similar t o the method developed by 
IVarshalovich & Levshakovl lll993h in order to infer the cosmo- 
logical variability of the proton-electron mass ratio fi from the 
analysis of molecular hydr ogen absorption lines. Argueing by 
analogy llLevshakovll2003l) . in the regime \Aa/a\ «: 1 we ob- 
tain the linear approximation 

Zi^Za+KaQi, (4) 

where z, - Aob^/Aoj and Qi = qi/dJoj denote, respectively, the 
observed redshift and the sensitivity coefficient corresponding 
to the line /, and the slope parameter Ka is given by 

K,^-2il+Za){Aa/a). (5) 

If Aa/a is non-zero, z, and Q, will be correlated and we will be 
able to estimate the slope and the intercept z„ from the lin- 
ear regression analysis of the position of the line centroids in an 
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absorption component. The accuracy of the regression analysis 
will be improved, if several absorption line samples are com- 
bined. In this case, the regression procedure can be generalized 
appropriately: 

^, = (Aff/a)(G,-G), (6) 

where Q denotes the mean sensitivity coefficient of the sample, 
and - (f -z,)/[2(l +z)] is the normalized redshift while z de- 
notes the mean redshift of an absorption component (in radial 
velocity space ^ = (S - m,)/[2(c + u)]). 

In order to determine the central position of several selected 
Fe II lines (see Sect. 14. 2> we basically follow the same strategy 
as described in Sect. 13.11 The only differences are that we do 
not incorporate Eq. Q and do not confine the radial veloci- 
ties to be the same for the lines in the selected components. 
We point out explicitly that even though we apply a parametric 
profile decomposition technique to determine the position of 
the line centroids, the RMM analysis is not tied to any specific 
modelling technique. In principle, the position of the line cen- 
troids can even be determined without doing any modelling at 
all (cf. lEevshakov 200J). Throughout the analysis we use the 
atomic data listed in Table [T] 

4. Results and discussion 

4.1. Standard MM analysis 

The optimized values of the model parameters and the stan- 
dard deviations provided by the covariance matrix of the ES 
are listed in Table |2l The decomposition of the Pen absorp- 
tion complex is quite evident (see Figs. ^ and IZ)- This con- 
trasts with the QSO absorption systems considered in former 
MM studies where many components ar e typic ally unresolved 
and many line profiles are saturated (see lMWFi Figs. 3 and 4). 
Clearly, the components with the most accurately defined line 
centroids (No. 1,12, 20-24, 26) are the most important in the 
MM analysis. Expectedly, the adequate profile decomposition 
is reflected in the formal accuracy of the result: {Aa/a}MM = 
(0.1 + 1.7) ■ 10"^ is the most stringent formal bound hitherto 
infered from an individual QSO absorption system. T he best 
formal accuracy achieved in former MM analyses (see lMwR 
Table 3) is exceeded by a factor of about three. 



Table 2. Standard MM analysis: optimized values and formal 
standard deviations of the radial velocity u, the line broadening 
velocity b, and the column density of the 29 components 
constituting the Fe ii absorption complex. The optimized value 
of {Aa/a} is {Aa/a)MM = (0.1 + 1.7) ■ 10"^. The zero point of 
the radial velocity corresponds to the redshift z - 1.1508 



No. w(kms ') fc(kms ') logA'(cm ') 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 



determination is - 0.005, i.e. only half a percent of the vari- 
ation among the relative displacement of lines is accounted for 
by the difference in sensitivity. Consulting the t statistic, the 
observed data support the null hypothesis (Aa/a) = at a sig- 
nificance level of 91 percent, whereas the support for the MEW 
result (Aa/a) = -5.7 ■ 10"^ is 12 percent. 



-534.35 + 0.07 2.07 

-530.62 ± 0.44 7.26 

-510.27 ±0.18 7.92 

-494.49 ±0.16 4.86 

-471.69 ±0.12 2.99 

-448.58 ±0.38 2.21 

-416.88 ±0.08 1.76 

-403.01 ±0.16 2.19 

-393.52 ±0.10 3.39 

-375.95 ±0.16 4.38 

-321.19 ±0.30 6.98 

-239.01 ±0.08 1.66 

-235.28 ±0.16 5.15 

-182.74 ±0.45 4.02 

-175.12 ±0.47 4.92 

-137.26±0.10 3.11 

-123.47 ±0.14 4.47 

-53.76 ±0.38 6.04 

-34.40 ±0.21 5.81 

-14.11 ±0.09 5.09 

-1.07 ±0.05 3.76 

8.71 ±0.06 3.46 

22.12 ±0.04 4.53 

37.33 ± 0.08 6.04 

59.75 ±0.16 11.3 

72.71 ±0.06 5.69 

88.98 ±0.24 9.16 

110.09 ±1.44 11.4 

134.09 ±0.43 4.08 



± 0.16 


12.29 ± 


0.02 


± 0.41 


12.06 ± 


0.03 


± 0.36 


12.14 ± 


0.02 


± 0.32 


1 1.91 ± 


0.02 


± 0.20 


1 1.91 ± 


0.02 


± 0.87 


11.15 ± 


0.06 


± 0.20 


1 1.91 ± 


0.02 


± 0.37 


1 1.79 ± 


0.02 


± 0.22 


12.12 ± 


0.01 


± 0.22 


1 1.77 ± 


0.02 


± 0.58 


1 1.52 ± 


0.04 


± 0. 1 1 


12.55 ± 


0.02 


± 0.15 


12.51 ± 


0.02 


± 0.42 


12.07 ± 


0.07 


± 0.40 


12.18 ± 


0.05 


± 0.18 


12.05 ± 


0.01 


±0.20 


12.00 ± 


0.02 


±0.59 


11.43 ± 


0.05 


±0.36 


11.77 ± 


0.03 


±0.10 


12.85 ± 


0.01 


±0.06 


13.56 ± 


0.01 


±0.10 


13.28 ± 


0.01 


±0.06 


13.54 ± 


0.01 


±0.14 


13.13 ± 


0.01 


±0.2 


13.48 ± 


0.01 


±0.08 


13.50 ± 


0.01 


±0.52 


12.51 ± 


0.02 


± 1.7 


11.70± 


0.08 


± 1.00 


11.19± 


0.07 



4.2. Regression MM analysis 

The optimized central positions of the lines considered in the 
RMM analysis are listed in Table |3] We do not consider com- 
ponent No. 22 since the central positions of the lines in this 
component are strongly correlated with those in No. 21 and 
we ignore component No. 24, because the central position of 
the Fell /11608 line is not accurately defined. The regression 
line (Fig. |5} indicates no correlation between the relative dis- 
placement of the lines, and the sensitivity coefficient Q: 
{Aa/a}RMM - (-0.4+1.9+2.7sys)-10"'', where the statistical er- 
rors of ( and the estimated errors of Q as well as the systematic 
uncertainties inherent to the wavelength calibration are propa- 
gated by means of Monte Carlo simulation. The coefficient of 



4.3. Systematic effects 

Contributing to the discu ssion o f potential systematic effects 
presented in ImWFI and ImEWI we recall attention to an- 
other important source of erro r. As formerly il l ustrate d by 
iLevshakov & D'Odoricol (Il995h and iLevshakovl ( 1 1994 . the 
presence of unresolved narrow lines with different optical 
depths can strongly affect the position of the line centroids in 
an ensemble of lines and result in a biased expectation value 
(Aa/a). This effect will be most noticeable in the case of op- 
tically thick lines. In fact, the expectation value {Aa/a}MM in- 
creases by 3.6 ■ 10"^ if we do not ignore the ensemble of satu- 
rated Fe II /12383 and /12600 profiles in the optimization, and the 
systematic increase of scatter in the normalized residuals with 
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Table 3. Regression MM analysis; optimized centroid positions (A) of the Feii lines between /11608 and /i2600 in the selected 
absorption components. Formal standard deviations are indicated in parentheses 



No. 


1608 


2344 


2374 


2383 


2587 


2600 


1 


3453.2908 (57) 


5032.9469 (24) 




5115.7159(17) 


5553.4531 (29) 


5582.4812(19) 


12 


3456.6993 (39) 


5037.9119(23) 


5102.9147 (36) 


5120.7677 (21) 


5558.9313(28) 


5587.9914(21) 


20 


3459.2892 (37) 


5041.6969 (25) 


5106.7434 (36) 


5124.6089 (22) 


5563.1057 (30) 


5592.1855 (24) 


21 


3459.4431 (11) 


5041.9145(12) 


5106.9734(12) 




5563.3452(11) 




23 


3459.7131 (13) 


5042.3034(10) 


5107.3695(16) 




5563.7764(11) 




26 


3460.2938(17) 


5043.1547(15) 


5108.2272 (21) 




5564.7129(16) 





4 - 




-4 



-0.04 -0.03 -0.02 -0.01 0.01 
Sensitivity Coefficient, Q-(3 

Fig. 3. Regression MM analysis of the lines listed in Table |3] 
Each component is represented by a different symbol. The re- 
gression line (Aa/Qr)RMM = -0.4 ■ 10"^ and its 68 and 95 
percent confidence limits are marked by the solid and dashed 
curves, respectively, whereas the dotted line indicates the MFW, 
result (Aa/a) - -5.7 ■ 10"*'. The median total errors of ( and 
Q are illustrated. 

increasing optical depth (see Figs.0]and|5lprovided in the on- 
line material) may be explained by the presence of unresolved 
narrow lines. 

5. Conclusions 

Our results strongly support the null hypothesis of a non- 
varying fine-structure constant, but do not contradict the MFW 
result (Aa/a) - -5.7 ■ 10"* at a significance level higher than 
88 percent. Nevertheless, we conclude: (i) The MM technique 
has the capability to provide stringent bounds on (Aa/a) even 
if the analysis is restricted to the Fe ii fines only, (ii) The RMM 
technique is illustrative and methodically more transparent than 
the standard MM technique. In addition, the regression analysis 
facilitates the consideration of systematic errors inherent to the 
wavelength calibration of QSO spectra, (iii) The accuracy of 
the individual Aa/a assessment is principally limited by sys- 
tematic errors inherent to the wavelength calibration. The ac- 
curacy attainable with high-quality QSO spectra recorded with 
the VLT/UVES is limited to 10"^. (iv) Optically thick profiles 
are susceptible to systematic effects biasing the expectation 



value (Aa/a). (v) The analysis of an extensive homogenous 
sample of Fe n absorption lines is inevitable and will provide 
an independent and crucial test of the MFW result. 

In particular, HE 0515-4414 is the brightest known inter- 
mediate redshift QSO in the sky and is therefore predestinated 
for spectroscopy with the new High Accuracy Radial velocity 
Planet Searcher (HARPS) operated at the ESO La Silla 3.6 m 
telescope. This very high-resolution spectrograph is specified 
to provide an eflicient wavelength calibration facilitating the 
performance of radial velocity measurem ents with an accuracy 
of better than 1 ms"' jPepe et al.ll20 02'). Presumably, HARPS 
will improve the accuracy of the individual Aa/a assessment 
by more than an order of magnitude. 
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Fig. 4. Normalized residuals of the optimized RMM model. For the purpose of residual analysis the calculation of the sum of 
residual squares takes into account the optically thick parts of the line ensemble. The systematic increase of scatter with increasing 
optical depth may be explained by the presence of unresolved narrow lines or by deviations from the Doppler profile resulting 
from e.g. kinematic effects due to rotational motion JProchaska & Wolfell997i) 
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Fig. 5. Number distribution of the normalized residuals rendered in Fig.@] The distribution is approximately normal and does not 
indicate any missing components 



